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•r:^ ' Abstract 

o: 

\^ . In view of the evaluation of the physics goals of a large Liquid Argon TPC, evolving 

^^ I from the ICARUS technology, we have studied the possibility of performing precision 

,^ ' measurements on atmospheric neutrinos. For this purpose we have improved existing 

Monte Carlo neutrino event generators based on FLUKA and NUX by including the 



Ph: 



Qh! 3-flavor oscillation formalism and the numerical treatment of Earth matter effects. By 






means of these tools we have studied the sensitivity in the measurement of 623 through 
the accurate measurement of fg's. The updated values for Attics from Super-Kamiokande 
and the mixing parameters as obtained by solar and KamLand experiments have been 
rS ' used as reference input, while different values of ^13 have been considered. An exposure 

C^ ' larger than 500 kton yr seems necessary in order to achieve a significant result, provided 

that the present knowledge of systematic uncertainties is largely improved. 

1 Introduction 

In the framework of the design study for a very large Liquid Argon TPC evolving from 
the present status of the ICARUS project^, it is important to start the discussion of the 
different physics goals that should be addressed by a multipurpose detector of this kind. 
The topic of atmospheric neutrinos is certainly one of the relevant chapters. Indeed, after 
the results coming mostly from Super-Kamiokande j2j, there is still a scientific interest in 
continuing the study of atmospheric neutrinos. First of all it is necessary to confirm the 
results from SK with a technology capable of a large reduction of experimental systematics 
with respect to water Cerenkov. In second place, many authors have pointed out the 
possibility of exploring subleading contributions in the oscillation matrix. This would 
allow, in particular, a possible precision measurement of 623 and the discrimination of 
the normal vs. the inverted hierarchy in the neutrino mass scale. These are tiny effect, 
and a very large Liquid Argon TPC with a sensitive mass of the order of several tens 
of kilotons seems necessary to reach the exposure necessary for these proposed analyses. 



The preliminary work here presented aims to investigate quantitatively the attainable 
performance of such a large LAr detector in the field of atmospheric neutrinos 

In section 1^ we describe in more detail the motivations of physics for this study and 
provide a general discussion of the effects we are going to discuss. In sectional we present 
the model implemented in our simulation to include three flavor oscillations and matter 
effects. Then, in section 0] we discuss the possible results with atmospheric neutrinos in 
the limit of a very high exposure. 

2 Atmospheric neutrinos and 3— F oscillations 

A large number of results regarding atmospheric [21 , solar |Hj, reactor [IJ [3] and accelerator 
neutrinos can be accounted for assuming that the ordinary 3 neutrinos have mass and mix 
among them jBj. In this framework, flavor oscillations are fully described by assigning 3 
mixing angles and 1 CP violating phase. We do not know the type of mass hierarchy, 
and 3 of these quantities are still unmeasured: The size of ^13, the deviation of ^23 from 
maximal mixing, and the CP violating phases 6 (see sect. |21for their definitions). The 
first unknown quantity is particularly important, for ^13 controls the size of three flavor 
oscillations. Although there are little doubts that such a Sv framework is appropriate to 
describe these data, it is also fair to say that any experimental indication that we have 
is appropriately described by 2 flavor (2F) oscillations and 3 flavor (3F) effects have been 
not yet seen. 

We recall that there are several theoretical approaches to understand neutrino masses, 
and possibly to predict unknown quantities as ^13. The problem of fermion masses is 
known to be a difficult one, and caution is necessary in interpreting any prediction or 
theoretical expectation. However, several models suggest that ^13 is not undetectably 
small (see e.g., |2j for SO(IO) models); rather, ^13 often happens to be close to the 
experimental limit of about 10°. Similarly, there are several models where deviations 
of ^23 from maximal mixing amount to several degrees. Presumably, these expectations 
can be considered in agreement with the common sense: in absence of any strong reasons 
to assume the contrary, ^13 and ^23 — 45° should be not too small. Certainly, experimental 
investigations of ^13 and the search of 3F effects in oscillations are among the most 
important goals still to be achieved and will have a profound impact on what we know 
about neutrinos. 

One of the most important results emerging from the 2F analysis of atmospheric 
neutrino experiments is that the ^23 mixing angle, which is defined in the range 0-7r/2, 
is found to be compatible with the full mixing value of 7r/4. The measurement of this 
mixing angle comes from the up-down asymmetry of muon-like events, and in particular 
of the higher energy events, as the Multi-GeV selection in Super-Kamiokande, where 
the direction correlation between neutrino and lepton is strong. Actually, such up/down 
ratio, which in the average results very close to 1/2, gives a measurement of sin^2^23 and 
therefore it is not possible in this way to distinguish the octant of 623- The present result 
is sin^2923 = 1 ± 0.01, which corresponds to siin?923 = 0.5 ± 0.05 In order to achieve a 
better understanding of the properties of neutrino mixing, it is of the utmost importance 
to measure precisely the value of ^23 and to measure its deviation from the full mixing 
value. 

The possibility of performing a different measurement emerges in the framework of 
the 3-flavor oscillation scenario. Beside the main part due to Vfj_ — Vr oscillations certainly 
we have further small effects of Ve oscillations with "solar frequency" and possibly, other 



fe oscillations with "atmospheric frequency" due to ^13. On top of these two effects, 
essentially of 2F character, there are also genuine 3F effects due to the interference between 
the various amplitudes of transition. 

The following analytical considerations help to understand the underlying physics. In 
general, one can write the formal solution for the amplitude Au of neutrino propagation 

as 

// Oil ai2 ai3 \ 

dtn^it)] = i?23A • 021 022 a23 A*i?*23 (1) 

\ 031 a32 033 / 

where i?23 are rotations by the angle ^23 in the 23 plane. The elements aij form an unitary 
matrix that depends on ^13 and 612 but not on ^23 or on 5 (of course, it depends also on 
Attic's and on the 'kinematical variables' E, L and 9z)- From the explicit form of the 
hamiltonian TCi, one easily shows that in the limit ^13 = 0, one has 013 = 023 = 0, so 
the effects due to 6*12 are described by 012 ('0i2-driven' or 'solar frequency' oscillations). 
Similarly, in the limit 612 = 0, one has ai2 = 023 = 0, so the effects due to ^13 are described 
by ai3 ('^13-driven' or 'atmospheric frequency' oscillations). 

In general, one can get aij from a numerical calculation. In the case of vacuum 
oscillations, however, the equations can be integrated and give for instance: 

ai2 = si2Ci2Ci3(e""^ - 1), ip = Amj^L/lE 

ai3 = si3Ci3[(e-**-l)-s?2(e-*''-l)], <^ = Aml,L/2E ^^> 

Here we are using the standard notation Sjj, cij for sinOij and cosOij respectively. It is 
evident that these amplitudes have properties just discussed in the limiting cases 612 -^ 
and 6*13 -^ 0. 

From eq. (^ one obtains expressions for the probabilities of survival or of conversion, 
e.g.: 

Pe^e = K-^!/)lll = 1 — |fll2l ~ 1^13! (o\ 

P^,^e = |(X)l2l = |C23 ai2 c'^ + S23 013^ ^ ' 

The interference term in the latter expression is a prime example of a genuine 3 flavor 
(3F) effect, that disappears in the limits 612 -^ and ^13 -^ 0. It should be noted that 
the dependence of this interference term on 013 is linear, thus it depends mildly on ^13. 
(It should be recalled that the CHOOZ or reactor limit on ^13, instead, scales with Of^.) 
Using eqs. (jHI and following the discussion in ^8^, let us introduce the flavor ratio 



r 



Fjl/F^. The variation of the electron neutrino flux F^, = Pe^eE^ + P^^eE?^ reads: 



4 - 1 = (rch - l)l«?2l + (^^23 - l)l«?3l + 2rs23C23Re[al^aue''] (4) 

e 

The first two contribution are small for low energy neutrinos, since r ~ 2 and the ^23 
mixing is close to maximal: C23 ~ S23 ~ l/v2- It should be noted that these two terms 
correspond largely to 2 flavor effects. In fact, from what told above we know that in the 
limit 6*13 = 0, ai3 = thus from the flrst of eqs. ((SJ |a^2l ~ ^ ~ -fe^e(Am^2' ^12) whereas 
in the limit 612 = 0, ai2 = thus \a1^\ = 1 — Pe-^e(A?7if3, ^13). The third term is the 
only one that it is not suppressed by the flavor ratio r ~ 2. In fact, this is exactly the 
interference term, namely, the genuine 3F effect. Note finally that the three terms work 
in different ways: the first one increases the electron flux if ^23 < 45°, the second one if 
^23 > 45°, the last, can increase or decrease depending on 6. In other words, the "solar" 
sector of the neutrino mixing (described by the AvTifg and 612 parameters) influence the 



rate of fe events (in particular in the SubGeV range) with respect to the no-oscillation 
case, with a different sign according to the value of ^23) even in case ^13 is very small or 
null. The effect vanishes totally instead for the maximum mixing case ^23 = vr/4. Since 
the solar neutrino and KamLand[2, 3] experiments have obtained a remarkably precise 
determination of two parameters ATTifg and 612, in principle the proposed measurement 
is possible, provided that systematics uncertainties are kept under control. 

The MSW effect [^ ("matter effects") is relevant in specific energy regions, and in 
particular this happens where vacuum term is comparable with the matter term. For 
instance, the vacuum term at Am^ = 10~^ eV^ and E = 0.5 GeV is Am? /2E w 5 • 10~^ 
eV^/MeV while the matter term at pe = 2.5 e" moles/cc is VlGpNe = 9 • lO""^ eV^/MeV 
(as usual, A^e = ^APe)- An interesting consequence is that, due to the MSW effect, ^13 
oscillations (driven by Am^^^ ~ 2.5 • 10~^ eV^) are amplified for £^ ~ 2 — 6 GeV. This 
could be of particular interest in long-baseline experiments, while, for statistical reasons 
connected to the steeply falling energy spectrum, atmospheric neutrino experiments are 
less sensitive in this energy range. Therefore in this work we focus on low energy events 
and neglect measurements aimed to measure ^13 and more in general MultiGeV events. 
We remark here that one of the most appealing advantage offered by the Liquid Argon 
technology with respect to water Gerenkov detectors, is the much larger efficiency in 
SubGeV investigation, in conjunction to the partial possibility of a better reconstruction 
of their kinematics. In fact an ICARUS-like detector can reconstruct recoiling protons 
down to a rather low energy threshold. 

There is another important measurement connected with 3F oscillations including 
matter effects which can be in principle performed using atmospheric neutrinos. This 
is the discrimination of normal from inverted mass hierarchy (the sign of Am\^). Large 
matter effects occur for neutrino if Amfg > and for antineutrinos if AvTT-fg < and 
this is again more significant in the MultiGeV region. As suggested in [lOj . a detector 
capable of distinguishing on an event by event basis v^ from V^ in the range from about 
3 to 10 GeV is ideal for this purpose. Such a separation capability can be performed 
by means of a magnetized detector: there is not yet a convincing design of large Liquid 
Argon coupled to a magnet. Furthermore, a low density detector has a low containment 
efficiency for MultiGeV muons. Therefore we also will not further discuss this suggestion, 
at least in this preliminary work which concerns the possibilities of further investigation 
of atmospheric neutrinos with a large Liquid Argon detector. 

3 The simulation code for 3F oscillations with 
matter effects 

The discussion of the previous paragraph introduces the need of calculating predictions 
taking into account with the maximum possible accuracy 3F oscillation coupled with 
matter effects. A numerical code has been designed for this purpose and we recall here 
the basic formalism that has been implemented for this purpose. 

The propagation of neutrinos or antineutrinos in a medium with electron density N^ 
that varies along the neutrino trajectory (parameterized by the time t) is described by 
Schrodinger-like evolution equations J) : 

Huit) = U ■ diag(m2) • C/t/(2S) + y^GpNeit) diag(l, 0, 0) 
njjit) = U* ■ diag(m2) . uy{2E) - V2GFNe{t) diag(l, 0, 
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Figure 1: Density profile of tlie Eartli according to tlie PREM model jTTj. 



where the neutrino mixing matrix includes the already mentioned three rotations through 
the 23, 13 and 12 planes with angles ^23) ^13 and 612 and a CP violating phase 5 



U = R23ARi3A*Ri2, A = diag(l, 1, e^' 



(6) 



In the following, we will discuss the case of neutrinos, however it is clear that the same 
discussion applies to antineutrinos inverting the sign of N^ and of the CP-violating phase 5. 



We use the notation Amf- 
whereas Aml^ 



7n^ — mf, and adopt the convention that Atti^j = ^''^soi ^ ^' 
dzAm^j^ if the mass hierarchy is normal or inverted respectively. In 
other terms, we declare that vs is respectively the heaviest or the lightest neutrino. 

In the code, the evolution equations above are solved resorting to numerical methods. 
Being concerned with atmospheric events, we consider neutrinos of energy E that reach a 
detector with zenith angle 6z after a path length L; if 9z > 90°, these neutrino traverse 
the Earth and therefore meet a non-zero electron density N^. We describe the Earth 
electron density by adopting the PREM model ^5. which gives the radial profile of the 
Earth's density as shown in Fig. ^ 

We solve the resulting evolution equations through a straightforward numerical 
integration.^ Given a neutrino (or antineutrino) , we calculate the probabilities that it 
will be detected as I'e, ^^ or Vj- and finally, sort the detected flavor according to these 
probabilities. 

This code has been coupled to a neutrino interaction generator (for both CC and 
NC) which makes use of the FLUKA I3. and NUX[T1] libraries. Quasi-elastic (QE) and 
deep inelastic scattering (DIS) are considered. The QE part is managed by FLUKA and 
according to the physics choice of NUX, the resonance region is considered in average. 



^We use the ODEINT subroutine given in Numerical Recipes [T^. with the yscal(i) variable set to unity. We 
obtained a numerical precision in the determination of the oscillation probabilities much better than one part 
per million for all neutrinos we propagated. This tight requirement is much more than what we actually need, 
but the required computational time is small when compared with the time to simulate an event. 



according to the duality principle, by the DIS part. FLUKA provides the nuclear 
environment for all reactions. Quantum and nuclear effects, including reinteractions in the 
nucleus, are simulated in detail. This event generator was originally developed within the 
ICANOE proposal!^ with the possibility of invoking 2- flavor oscillations (full mixing) 
in vacuum. Now the present version has been upgraded so to include the general 3- 
flavor oscillation code coupled to the above described numerical solution of the differential 
equation for the transport of neutrino amplitude through the Earth. 

Oscillation parameters are provided by means of a data card file {osc. datacards), where 
Atti^ values must be given in eV'^ and mixing angles in degrees. The CP violating phase 
(5 = (0 — 360°) is also considered, and the user can choose the "direct" or "inverse" 
hierarchy. Matter effects can be switched off. A typical example of the osc. datacard file 
is reported below: 

#Data cards for atm. neutrino oscillation code 

# Thetal2 in degrees 
THETA12 32.3 

# Thetal3 in degrees 
THETA13 0. 

# Theta23 in degrees 
THETA23 45 . 

# DeltaM12**2 in eV**2 
DMSQ12 8.2D-5 

# DeltaM23**2 in eV**2 
DMSQ23 2.6D-3 

# CP violation phase in degrees 
DELTA . 

# Hierarchy: normal: 1, inverted: -1 
HIERARCHY 1. 

# Activate matter effects: on: 1, off = 
MATTER 1. 

# end of data cards 
END 

These drivers require as input the neutrino fluxes. For atmospheric neutrinos we start 
from the "FLUKA" fluxes of 2001 J16j. which are calculated for 3 specific geographic 
locations. For the purpose of this work, we limited ourself to the LNGS case. 

As a default choice the interactions are performed against argon nuclei, but any other 
nucleus can be considered. The natural isotopic composition of any given nucleus is 
automatically provided by the FLUKA database. These generators produce in output 
the detailed kinematical information of the initial, intermediate and final state of the 
interaction and, optionally, the binary input for a full simulation of detector response as 
deriving from ICARUS. 

4 Atmospheric neutrinos in a large Liquid Argon 
TPC 

We have simulated a single experiment with an high exposure, simulating the statistical 
sample of CC interactions of atmospheric neutrinos at Gran Sasso expected in 1000 



kton yr. We have considered, beyond the unoscihated case, three cases for ^23 well inside 
the present limit from atmospheric experiments: 40°, 45° and 50°. The fit results from 
the review of ref.|BI have been used as reference for A?tt-23, "^"^12 ^^^ ^12 (2-5 10~^ eV^, 

8.0 10~^ eV^, 34°), while ^13 has been varied between 0° and 10°, {i.e. about the existing 
limit from the CHOOZ experiment ;5 J. For the work described here, the CP violation 
phase 6 has been left to 0° and only the direct hierarchy of masses has been considered. 

4.1 Distribution in zenith angle and L/E 

Before entering into the details of the measurement of the octant of ^23 it is worthwhile 
considering the "standard" measurements with atmospheric neutrinos. For instance, FigEl 
shows the expectation for the zenith angular distribution of muon neutrino events, when 
only the lepton direction is measured, for SubGeV and MultiGeV regions in the case of 
full mixing (^23 = 45°). The two extreme values of ^13 are considered in the oscillation 
predictions. 

This plot already gives a quantitative idea of the difficulty of measuring ^13 with 
atmospheric neutrinos. However, a Liquid Argon TPC gives the chance of performing 
an improved measurement of the neutrino direction. This is already possible, in part, 
in the SubGeV range. Notwithstanding the intrinsic smearing due to Fermi motion, a 
better correlation with original neutrino direction can be maintained, as shown in Fig|3J 
where the kinematic information of possible recoiling proton has been added, assuming 
a detection threshold of 50 MeV for the kinetic energy of an hadron. The expected 
enhancement around the horizontal direction is clearly visible. Despite the reduction in 
statistics, the separation between the two extreme cases of ^13 is slightly improved. 

The above considerations suggest that the results from a large Liquid Argon TPC can 
be used to perform a significant identification of the oscillation pattern by means of the 
plot of the ratio of event rate with respect to the non oscillated prediction as a function 
of log L/E. This is reported in Fig^ (L is measured in km, while E is in GeV). There we 
have followed the usual method of mirroring the direction of downward going events to 
obtain a model independent prediction for upward going neutrinos. In the upper panel 
of Fig^ only the lepton direction has been considered, while in the bottom panel the 
kinematics of the hadronic system (for particles above threshold) has been included: the 
significance in pattern recognition is clearly improved. 

We estimate that, for this measurement, an exposure of 500 kton yr would be already 
sufficient to provide a very significant result. 

4.2 Low energy electron events and 623 

We enter now the topic of the precision measurement of ^23- According to the discussion of 
section|2l we expect a variation of the rate of electron neutrino events (and in particular of 
SubGeV electron neutrinos) as a function of 623. In the simple case of vacuum oscillations 
and ^13 we expect: ANg oc (1 — N^/N^siv?92'i), where N^ and N^ are the unoscillated 
/x-like and e-like event rates. The situation is slightly more complex for ^13 = and when 
matter effects are present. Our complete simulation gives the behavior shown in Fig. [Sj 
where the deviation from no-oscillation prediction of event rate of SubGeV e-like events 
(in events/kton yr) is plotted vs siv?923. 

As expected, the relative differences in rate are very small and a very large exposure 
is needed. The absolute difference in rate of SubGeV e-like events is shown in Fig. El for 
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Figure 2: Zenith angle distribution for SubGeV (top) and MultiGeV (bottom) /i-like events. 
The non oscillated and oscillated predictions for 6*23=40° and two different values of 6*13 are 
reported. Here only the lepton direction has been used. 
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Figure 3: Zenith angle distribution for SubGeV /i-like events where both lepton and recoihng 
proton kinematics has been reconstructed (Efcj„ > 50 MeV). The non oscillated and oscillated 
predictions for ^23=40° and two different values of Oi^ are reported. 



^13 = and for two values of O23 syininetric with respect to 45° (only the lepton direction 
is considered here). 

If we take the ratio with respect to the non oscillated prediction, we get the result 
plotted in Fig|7| 

If we conservatively fit the SubGeV distribution with a constant (namely this is 
equivalent to integrating the angular distribution over the whole range of the zenith angle) 
the difference in the ratio N^/N^ between 623 = 40° and 50° is at the level of 0.037 ± 
0.006. For relatively large values of ^13, the significance in the difference becomes smaller, 
as shown in Fig|Hl 

Of course it is hard to believe that one could rely on the absolute level of A'^e prediction: 
flux normalization remains one of the most important uncertainties. A better analysis is 
obtained by the double ratio: {Nf./N^)/{N^/Nf^). The variation of /i-like events has the 
opposite trend of A'^e, and in first approximation all common uncertainties contributing 
both to A'^e and A'^^ cancel out. Here the most important uncertainties concern the 
normalization of the primary cosmic ray flux and the knowledge of i/-Nucleus cross 
sections. In good approximation these affect e and /i flavor almost at the same level 
(there are mass dependent terms which are of course non negligible at very low energy). 
As an example of the possible expectation, we get the behavior as a function of the zenith 
angle which shown in Fig|ni(when ^13 = 0°). 

Now the important topic, as far as systematic uncertainties are concerned, is the 
accuracy on the prediction of the N^/N^ ratio. This topic has been already debated in 
the context of a dedicated workshop stimulated by Suoer-KamiokandefTT] . One of the 
reported conclusions was that the measurement of 623 octant can be done achieving an 
exposure of at least 20 years of runs of Super-Kamiokande in order to distinguish (at the 
level of Ax^ ~2) between the two mirror values corresponding to sin 2^23=0.96 (~39°-^ 
51°), with the present level of systematics. An effort to improve the accuracy on the 
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Figure 4: Ratio of event rate with respect to the non oscillated prediction as a function 
of log L/E. L is measured in km, while E is in GeV. Top panel: only lepton direction is 
reconstructed. Bottom panel: both lepton and charged hadrons are reconstructed (Efcj„ > 50 
MeV). 




Figure 5: Deviation of the SubGeV e-like event rate (in events/kton yr) from the no-oscillation 
prediction as a function of siv?62z. 
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Figure 6: Absolute prediction of the zenith angular distribution of SubGeV e-like events for 
two different values of 6*23 when 6*13 = 0°. 
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Figure 7: N^./N'^ ratio as a function of the zenith angle of SubGeV e-hke events for two different 
values of ^23 when 6'i3 = 0°. 



prediction on A^^/A^g at a level better than 1% would be very helpful. 



5 Summary and discussion 



A very large LAr TPC allows to detect low energy neutrinos (above ~ 50 MeV) with null or 
negligible experimental systematic error. Therefore, in principle, such a detector can give 
new important contributions to neutrino physics, also by means of a new investigation of 
atmospheric neutrinos. An exposure of at least 500 kton yr is however necessary in order 
to give new significant contributions. The sector of SubGeV fg, which is particularly 
suitable for an ICARUS-like detector, offers in particular the possibility of performing 
new precision measurements. In this work we have performed a preliminary investigation 
about the possibility of the measuring how close is ^23 to the full mixing value of 45°. In 
principle this measurement is affordable if two conditions are met: an exposure larger or 
equal to 500 kton yr and a reduction of systematic uncertainties about neutrino fiuxes, 
where the most important item is the knowledge of the N^/N^ ratio. 
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Figure 8: Ne/N'^ ratio as a function of the zenith angle of SubGeV e-hke events for two different 
values of 6*23 when 61^ = 5° (top) and ^13 = 10°. 
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Figure 9: Double ratio {Ne/N^) / {N^/N^) as a function of the zenith angle of SubGeV e-like 
events for two different values of 6*23 when 6*13 = 0°. 
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